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r— I . Abstract 
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Due to the momentum transfer r = p — p' from the initial proton to the final, the "asym- 
metric" matrix element (p'\G . . . G\p) that appears in the pQCD description of hard diffrac- 
tive electroproduction does not coincide with that defining the gluon distribution function 
f g (x). I outline a lowest- twist pQCD formalism based on the concept of double distribution 
F g (x,y), which specifies the fractions xp, yr, (1 — y)r of the (lightlike) initial momentum 
p and the momentum transfer r, resp., carried by the gluons. I discuss one-loop evolution 
equation for the double distribution F g (x,y;/j 1 ) and obtain the solution of this equation in 
a simplified situation when the quark-gluon mixing effects are ignored. For r 2 = 0, the 
. ^ | momentum transfer r is proportional to p: r = (p, and it is convenient to parameterize the 

matrix element (p — r\G . . . G\p) by an asymmetric distribution function Ff{X) depending 
on the total fractions X = x + yC and X — £ = x — (1 — y)( of the initial proton momentum p 
carried by the gluons. I formulate evolution equations for (X), study some of their general 
properties and discuss the relationship between Ff{X), F g (x,y) and f g {x). 



1 Also Laboratory of Theoretical Physics, JINR, Dubna, Russian Federation 



1. Introduction. As shown in ref.p], at high virtualities of the virtual photon 7*, one can apply 
pQCD factorization to study the process of hard diffractive exclusive electroproduction of vector 
mesons 7* + p — > V + p' . In the approach of ref.|IJ, the non-perturbative information related to 
the proton is described by the matrix element of a two-gluon operator approximated by the gluon 
distribution function f g (x). However, as noted in ref.pl], due to the momentum transfer r = p — p' 
from the initial proton to the final, the two gluons carry, in fact, different fractions of the original 
proton momentum, i.e., the matrix element of the gluonic operator in this case does not coincide 
with that defining the gluon distribution function. As emphasized by X.Ji ||, one should deal in 
this case with a new type of functions (he calls them "off-forward parton distributions") which 
differ from f g (x), even if t = r 2 vanishes. My goal in this letter is to develop a modified pQCD 
approach for the hard diffractive electroproduction, which takes into account the effects related 
to the momentum transfer. In this formalism, the basic function describing the gluon content of 
the "asymmetric" matrix element (p — r\ . . . \p) is the double distribution F g (x,y), which specifies 
the fractions xp, yr, (1 — y)r of the initial proton momentum p and the momentum transfer r, 
resp., carried by the gluons|. With respect to x, the function F g (x,y) looks like a distribution 
function while with respect to y it behaves like a distribution amplitude. Since the logarithmic 
scaling violation is an important feature of the gluon distribution function in the low-x region, 
I discuss the evolution equation for the double distribution F g (x,y; fi). The relevant evolution 
kernel R gg (x, y; £, rf) produces the GLAPD evolution kernel P gg {x/^) [§], ||, || when integrated over 
y, while integrating R gg (x, y; £, rj) over x gives the expression coinciding with the evolution kernel 
Vgg{%),rj) for the gluon distribution amplitude. I construct the solution of the one-loop evolution 
equation for the double gluon distribution in a simplified situation when quark-gluon mixing effects 
are neglected. For t — and vanishing hadron masses, the momentum transfer r is proportional to 
p. r = (-p anc ^ f or reasorL) it is convenient to parameterize the matrix element (p — r\G . . . G\p) 
by the asymmetric distribution function (X) specifying the total fractions Xp, (X — Qp of the 
initial hadron momentum p carried by the gluons^. I formulate equations governing the evolution 
of the asymmetric distribution function T®(X) and discuss the relationship between this function, 
the double gluon distribution F g (x,y) and the usual gluon distribution function f g (x). 

2. Double distributions. The amplitude for the elastic electroproduction process 7*p — > p'V 
depends on the momentum p of the initial proton, the momentum transfer r = p — p' and the 
momentum q of the produced vector meson. We will consider the limit in which one can neglect 
the squares of the meson q 2 = my and proton p 2 = m 2 masses compared to the virtuality 
— Q 2 = (g — r) 2 of the initial photon and the energy invariant p ■ q. Thus, we set p 2 = and 
q 2 = 0, and use q and p as the basic Sudakov light-cone 4- vectors. In the diffractive region, the 
invariant momentum transfer r 2 = t is small, and we actually will analyze the limit t = 0. Note 
that in this case the on-shell condition p' 2 = (p — r) 2 = p 2 results in the requirement (p ■ r) = 

2 Originally, the double distributions for the asymmetric matrix elements of quark operators were introduced in 
ref. in application to virtual Compton scattering. 

3 The asymmetric distribution functions are similar to, but not identical with the t — * limit of the off-forward 
parton distributions introduced recently by X.Ji ||. 
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which can be satisfied only if the two lightlike momenta p and r are proportional to each other: 
r = Cp., where ( = Q 2 /2(p ■ q) is the Bjorken variable which obeys < £ < 1. 

The leading contribution in the large-Q 2 , fixed-C limit is given by the diagrams shown in Fig.|I| 
The long-distance dynamics is described there by the vector meson distribution amplitude fv{ T ) 
and the asymmetric matrix element of the light-cone gluonic operator 

(P " r | G« a ( Zl )E ab (z u z 2 - A)G b ua {z 2 ) \ p) (1) 

in which the 4- vectors Zi,z 2 specifying the location of the two gluon vertices are separated by 
lightlike intervals z\ = 0, z\ = from the virtual photon vertex. The factor E a b(zi, z 2 ; A) is the 
usual P-exponential of the gluonic A-field along the straight line connecting z\ and z 2 , and the 
indices a,b = 1, ... 8 denote the gluon color. 

Q- r xq q 



xp+yr 




<ZZ>> <_> -CLJ^ <ZS^ 

P p-r 



a) b) c) d) 

Figure 1: Diagrams contributing to hard diffractive electroproduction of vector mesons. 

The gluon momentum in this matrix element originates both from the initial hadron momentum 
p and from the momentum transfer r, so I write it as xp+yr parameterizing the asymmetric matrix 
element of the light-cone gluon operator by the double distribution F g (x,y) 

(p-r\ Zfl z u G; a (0)E ab (0, z; A)G b av {z) \ p) |, 2=0 (2) 

= u(p - r)zu{p) {z ■ p) £ jT 1 ^ (e-toM-M") + e «(p*)-*S(«)^ 0{x + y < 1) F g (x, y) dx dy. 

Here and in the following I adhere to the convention y = 1 — y, x — 1 — x, etc., for momentum 
fractions and use the notation 70,-2" = z. 

Due to the spectral properties x>0,y>0,x + y<l (this can be proved for any Feynman 
diagram using the approach of ref. ||), both the initial active gluon and the spectators carry 
positive fractions of the initial hadron momentum p: (x + (y) for the gluon and (x — (y) > 
y(l — (■) > for the spectators. On the other hand, the fraction of the p-momentum carried by 
another gluon is given by (x — y() and it may take both positive and negative values. 

The usual gluon distribution function xf g (x) corresponds to the limit r = 0. Hence, the double 
distribution F g (x,y) satisfies the reduction formula: 

rl—X 

/ F g (x,y)dy = xfg(x) . (3) 
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3. Asymmetric distribution functions. Since r = (p, the variable y appears in eq.(|2]) only in 
the combinations x + yC, = X and x — yC, = X — (, where X and (X — Q) are the total fractions of 
the initial hadron momentum p carried by the gluons (cf. 0). Introducing X as an independent 
variable, we can integrate the double distribution F(X — y(,y) over y to get 



min{X/C,X/C} 



F(X-y(,y)dy, 



(4) 



where ( = 1 — (. Since £ < 1 and x + y < 1, the variable X satisfies a natural constraint 
< X < 1. 




Qp 




(1-Y)r 



P-r 



a) 



b) 



Figure 2: Momentum flow corresponding to the asymmetric gluon distribution function in two 
regimes: a) X > ( and b) X = Y( < (. 



In the region X > ( (Fig. 2a), where the initial gluon momentum Xp is larger than the 
momentum transfer r = (p, the function JF^(X) can be treated as a generalization of the usual 
distribution function xf g (x) for the asymmetric case when the final hadron momentum p' differs 
by (p from the initial momentum p (but remains collinear to it). In this case, T®{X) describes 
a gluon going out of the hadron with a positive fraction Xp of the original hadron momentum 
and then coming back into the hadron with a changed (but still positive) fraction (X — ()p. The 
Bjorken ratio ( serves here as an external parameter specifying the momentum asymmetry of the 
matrix element. Hence, one deals now with a family of asymmetric distribution functions JF^(X) 
whose shape changes when ( is changed. The basic distinction between the double distributions 
F(x, y) and the asymmetric distribution functions F^(X) is that the former are universal functions 
which do not depend on the momentum asymmetry parameter (, while the latter are explicitly 
labelled by it. When ( — > 0, the limiting curve for J-^{X) reproduces the gluon distribution 
function: 

J^ =0 (X)=Xf g (X) . (5) 

Another region is X < £ (Fig. 26), in which the "returning" gluon has a negative fraction 
(X — C) of the light-cone momentum p. Hence, it is more appropriate to treat it as a gluon 
going out of the hadron and propagating together with the original one. Writing X as X = Y(, 
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we immediately obtain that the gluons carry now positive fractions Y(p = Yr and, respectively, 
(1 — Y)r = Yr of the momentum transfer r. Hence, in the region X = Y( < £, the asymmetric 
distribution function looks like a distribution amplitude ^ (V) for a two-gluon state with the total 
momentum r = (p: 

* f (y)= f Y F((Y-y)C,y)dy. (6) 

w 

The asymmetric distribution function can be also defined directly through the matrix element 

(p' | ^^(0)^(0, z; A)G b au {z) \ p)\ z , =Q (7) 
= u(p')zu(p) (z -P) f\ (e' zXipz) + e l{x ~^ z) ) T 9 C {X) dX. 

To re-obtain the relation between T^{X) and the double distribution function F g (x, y), one should 
combine this definition with eq.(0). 

4- Leading contribution. The parameterization (|7|) can be used as a starting point for con- 
structing a QCD parton-type formalism. The only problem is that our gauge-invariant definition 
of the gluon distribution is in terms of the field strength tensor G^ u while the usual Feynman rules 
involve the vector potential A^. A possible way out is to utilize the light-cone gauge q^A^z; q) = 
in which A^ can be expressed in terms of 

POO 

A^z- q) = q u G^{z + aq) da, (8) 
Jo 

so that the definition (|7|) can be applied directly. This gives 

ll\ A a( \ A a( \ I vi U(p')qu(p) f PnQv+PvlA / n x 

(p \A ll {z 1 ;q)A v {z 2 \q) |p)|^=o,^=o = 2( ; , \-~9^ + f . \ ) ^ 



x 



2(q ■ p) V (p-q) 

1 ( -iX(pzi)+i(X-C)(pz 2 ) , p »(A--0(p«i)-iX(p«) , \ ^Hj jy 

V 6 +e J X (X-C + ^e) 



In ref . [J] , the amplitude of hard diffractive electroproduction was calculated for the longitudinal 
polarization of both the virtual photon (e^* = (g M + (p^/Q) and produced meson (ey = q^/my). 
In this case, we obtain 

f 1 Mr) f 1 F 9 c (X)dX 



TLL(p,q,r)~l—^ ^dr ^ ' (10) 

Qm v Jo tt Jo A(A — Q + ie) 



where a/1 — £ comes from u(p') = a/I — ( u{p) and <£v( r ) is the distribution amplitude of the 
longitudinal vector meson. The ampitude has the imaginary part due to the factor 1/(X — ( + ie): 



71 QQmy s JO TT 



In ref.|l[], the gluonic matrix element was approximated by the gluon distribution function f g (C)- 
To get our result from that of ref.|j], one should substitute there f g (() by a/1 — C-^ (C)/C 
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Though the asymmetric distribution function T^{X) coincides with Xf g (X) in the limit ( = 0, 
these two functions differ in the general case when ( ^ 0. Furthermore, the imaginary part appears 
for X = (, i.e., in a highly asymmetric configuration in which the second gluon carries a vanishing 
fraction of the original hadron momentum. Hence, one cannot exclude the possibility that ^(C) 
visibly differs from the function (f g (() which corresponds to a symmetric configuration in which 
the final gluon has the momentum equal to that of the initial one. 

To get a feeling about the interrelationship between T^{X) and Xf g (X), let us consider a toy 
model F mod (x,y) = A(n + 1)(1 — x — y) n for the gluon double distribution. Then xf mod (x) = 
A(l -x) n+1 while 



X) n+1 - (1 - X/C) n+1 ) 6(X < C) + (1 - X) n+1 6(X > c) 



;i2) 



Hence, (f mod (() = A{1 - () n+1 while ^ od {() = A(l - C) n , i-e., the two functions are rather 
close to each other for small (. This model also reveals a characteristic feature of the asymmetric 
distribution function jF™ od (X): the parameter ( specifying the momentum asymmetry of the 
gluonic matrix element serves as a boundary between the two regions X < ( and X > ( in which 
J-^(X) is given by different analytic expressions. An important property of T^ lod {X) is that it 
rapidly varies in the region X< ( and vanishes for X = 0. Since X = can be arranged only 
when both x- and ^-parameters of the double distribution F(x, y) are set to zero, the fact that 
jF™ od (0) = is quite general. Note, however, that the limiting curve F™3(X) = (1 - X) n+1 does 
not vanish for X = 0, i.e., the limits ( — > and X — > do not commute. 

For this reason, if ( is small, the substitution of ^(X) by Xf g (X) may be a good approxi- 
mation for all X-values except for the region X< (. However, the imaginary part is given by the 
value of Fc,{X) at the point X = ( belonging to this region and it is not clear a priori how close 
are the functions ^(C) an d (fg{()- 

5. Evolution of the double distribution. On the light-cone, the matrix elements have ultraviolet 
divergences, which are removed by subtraction prescription characterized by a scale ji: F g (x, y) — > 
F g (x,y;fi). Under renormalization, the gluonic operator 

O g (uz, vz) = Zil z v G a m {uz)E ah {uz, vz; A)G h va {vz) (13) 

mixes with the flavor- singlet quark operator 

i _ _ 

O q (uz, vz) = - ^2(ip q (uz)zE(uz, vz; A)4> q {yz) - ip q (vz)zE(vz, uz; A)vj q {uz)) . (14) 

1 q 

For simplicity, we will ignore here the quark-gluon mixing and analyze below the evolution of 
the double gluon distribution F g (x,y; /i) corresponding to the "quenched approximation". In this 

case 

d ~ f 1 f 1 

v~djf^ x > mti } = J d tj R 9g( x 'y'^^;g^))Fg(^v^)dv- (15) 
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The easiest way to get explicit expressions for R ab (x,y;£,r];g) is to use the Balitsky-Braun 
evolution equation M for the light-cone operators^] 



fi—O a (0,z)= / / BgJu, v)O b {uz, vz) 6(u + v < 1) du dv . 
d\i Jo Jo ^ 

where a,b = g,Q and M 



(16) 



B, 



( u ,v) = — N c ( 4(1 + 3uv-u-v) + J^-6(u)5(v 
7r V 2iV c 



yy 



+S(u) 



v Jo z 



+ 6{v) 



U - - 5{u) t *Z 



u 



z 



(17) 



Here, (3q = \ \ — |iVy is the lowest coefficient of the QCD /3-function. 
Our kernel R gg (x, y\ £, 77; g) is related to the B gg (u, t>)-kernel by 



This gives 



a s 1 f 

Rgg(z,y,Z,ri;g) = — N c -{ 4[x/£ + 3(y - rjx/^)(y - WO] 0(0 < < min{?//77, y/fj}) 



7T ^ 

^(0<x/^< l)(a;/0 2 



(19) 



-6(x/£-y/rj) + -5(x/£-y/r]) 

T} 7] 



+ S(l-x/OS(y-v) 



2N r . 



1 dz 
z 



As usual, the divergent integral provides the regularization for the singularities of the kernel for 
x = £ (or y = rj). It is easy to verify that the kernel R gg (x, y; £, r\\ g) has the property that x+y < 1 
if £ + f] < 1. Using the expression for R gg (x, y; £, 77; g>) and the explicit form of the GLAPD kernel 
P gg (x/£) m |[, one can check the reduction formula 



l-x 1 

R 99 (x,y;^r];g)dy = -P gg (x/i). 



(20) 



Integrating R gg (x,y;£,,r); g) over x one should get the evolution kernel for the gluon distribution 
amplitude 



1-2/ 



Rggfa V, £, m 9)dx = V gg (y, rj- g). 



(21) 



To solve the evolution equation, we apply first the standard trick used to solve the 
GLAP equation: integrate x n R gg (x,y; £,77; g) over x. Using the property R gg (x,y;£,r);g) = 



R g f(y,r,-g)F^(r,^)dr,. 



(22) 



4 Instead of the original kernels K a b(u,v) from ref.|| we prefer to use the kernels B a b{u,v) = —K a b(u,v) which 
have the symmetry property B a b(u,v) = B a b(v,u) . 
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where F^ n \y; u) is the nth x-moment of F g (x, y; u) 



Ff\y\li) = fx n F g (x,y,fi)dx 



(23) 



and the kernel R g n g \y, 77; g) is given by 



n+2 



4 12 
n + 2 n + 1 



12 



+ 



1 



2N r 



1 d2 

z 



n + 2 

{y ^y,v 



(yrj + ny) 
v}\- 



12 



n + 3 



vn 



(24) 



It is straightforward to establish that R g T g > (y,r]; g) has the property 

R$(y, m g)w n {n) = R<$ (77, y- g)w n (y), 

where w n (y) = (yy) n+2 . Hence, the eigenfunctions of R g T g > (y,r]; g) are orthogonal with the weight 
w n (y) = (yy) n+2 , i.e., they are proportional to the Gegenbauer polynomials C£ +5 ^ 2 (y — y) (cf.fllOl, 



|TT| and refs. |T2| , |nj where the general algorithm was originally applied to the evolution of gluonic 
distribution amplitudes). Now, we can write the general solution of the evolution equation 



Ft\y, aO = (yy) n+2 E ^ k c n k + " 2 {y - y) [iog(/VA)P* 



(25) 



k=0 



where the anomalous dimensions are given by the eigenvalues of the kernel R^{y,n; g): 



mi 



(n) 
Ik 



2N r 



1 



1 



k+n+l j 

(fc + n)(fc + n + 1) (fc + n + 2)(/c + n + 3) + ^ j 



(26) 



They coincide with the standard anomalous dimensions ^y gg (N) |TJ], f5|: 7^ = 7 gg (n + + 1). 
Note, that 7° is formally given by the negative infinity, while all other anomalous dimensions are 
finite and non- negative. Hence, in the formal u — > 00 limit we have F n (y,n — > 00) = for all 
n > 1. This means that 

F fl (z,y;/x -> 00) ~ 5(x)(yy) 2 , 

i.e., in each of its variables the limiting function F g (x, y; u — ► 00) acquires the characteristic asymp- 
totic form dictated by the nature of the variable: S(x) is specific for the distribution functions 
Hl4] , |15| , while the (yy) 2 -foim is the asymptotic shape for the lowest-twist gluonic distribution 
amplitudes [12|, It is easy to see that if the double distribution has the asymptotic form 
F as (x,y) = l20C5(x)y 2 (l-y) 2 , then xf as (x) = C 5(x) while F% S {X) = 120 C X 2 (l - X/() 2 /(. 
Note that in this case .F° S (C) = 0, i.e., the function which determines the magnitude of the 
imaginary part of Tll vanishes. 
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6. Evolution equations for asymmetric distribution functions. Introducing the asymmetric 
distribution function JF^(X) for the flavor-singlet quark combination ([14]) 

d 



(p-r\O Q (0,z)\p)\ z 2 =Q = iu(p-r)zu(p) j (e~ lX ^ - e l{x ~^) (X) dX (27) 

and using eq.flIE|), we obtain a set of coupled evolution equations for !F 9 (X) and J^{X): 

^^(X; /i) = £ £ Wf(X, Z; g) f*{Z; /i) dZ , (28) 

where a and b denote g or Q. To relate the functions W£ b (X, Z; g) to the Balitsky-Braun evolution 
kernels || 

Bqq(u,v) = ^C F (l + 5(u)[v/v] + + 6(v)[u/u} + - ^(«)*(«)) , (29) 
(u,v) = — C F (2 + 5(u)5(v)) , B Q Ju,v) = —N f (l+4uv-u-v) (30) 

7T V / 7T 
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(B gg (u,v) was displayed earlier by eq.(|T^)), it is convenient to introduce first the auxiliary kernels 
Mf(X,Z;g): 

Mf(X, Z) = f 1 f 1 B ab (u, v) 6(X -uZ + v(Z - Q) 9{u + v < 1) du dv . (31) 
s Jo Jo 

In terms of these kernels, we have 

W% a (X, Z) = Mf (X, Z) , W^ Q (X, Z) = M C QQ (X, Z), (32) 

W 9 C Q {X, Z) = f 1 Mf(X, Z) dX , W? 9 (X, Z) = 4tf M? 9 (X, Z) . (33) 
Jx dX 

Note that the expressions for the kernels describing the quark-gluon mixing are slightly more 
involved than those for the diagonal ones. The reason is that the definition (|7|) of the gluon 
distribution has an extra (z ■ p) factor which generates a derivative acting on T 9 {X). 

Integrating the delta-function in eq . (|3l"D , one obtains four different types of the ^-functions, 
each of which corresponds to a specific evolution regime for the asymmetric distribution functions. 
In particular, for W 99 (X, Z\ g) we obtain 

W 99 (X, Z;g) = ^ jf* dv B gg (\l - X/Z - v(l - (/Z)} , (34) 
x{ 9(Z > X > C) (o < v < \~_f/z ) + °( Z ~ C - X) 0(0 - v ~ X/C) 



+0{X < C) 6{Z < C) 



6(x <z)6(o<v< x/o + 6(x > z) e I < « < x/c 



For the first two terms in this sum, the original fraction Z is bigger than the momentum 
asymmetry parameter (. In this region of Z's, the resulting fraction X cannot be increased by 
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the evolution: in both terms we have X < Z. Such a situation is typical for the evolution of 
distribution functions. 

The last two terms in eq.(|34]) correspond to the region where the original fraction Z is smaller 
than (. In this case, the evolution can either decrease the fraction (X < Z for the first term in 
the square brackets) or increase it ( X > Z in the second one). However, if the initial fraction Z 
is less than (, the evolved fraction X cannot be larger than ( or, what is the same, the parameter 
Y = X/( specifying the fraction Yr of the momentum transfer r carried by this gluon cannot 
exceed 1. This property is a characteristic feature of the evolution of distribution amplitudes. 

Qualitatively, the evolution of the asymmetric distribution functions proceeds in the following 
way. Due to the GLAP-type evolution, the momenta of the partons decrease, and distributions 
become peaked in the regions of smaller and smaller X. However, when the parton momentum 
degrades to values smaller than the momentum transfer r = (p, the further evolution is like that 
for a distribution amplitude: it tends to make the distribution symmetric with respect to the 
central point X = Q/2 of the (0, Q segment. In two extreme cases, when ( = or ( = 1, the 
evolution is more trivial. For £ = 0, ^(X) reduces to a usual distribution function governed by 
the GLAP evolution equation, while for ( = lwe always have Z < ( for the initial fraction and the 
function J-^(X) experiences a purely Brodsky-Lepage evolution. In other words, W^ (X, Z; g) = 
P ab (X/Z, g)/Z and W^(X, Z; g) = V ab (X, Z; g) Q. 

7. Conclusions. In this letter, I demonstrated that one can describe the asymmetric matrix 
element ((p\G . . -G\p) either by the universal double distribution F g (x,y) or by the asymmetric 
distribution function ^(X) which explicitly depends on the momentum asymmetry parameter 
( and specifies the total fractions X and X — ( of the original hadron momentum p carried by 
the gluons. Using ^(X) gives a formalism that looks very similar to the standard QCD parton 
approach, in which the gluon content of the hadron is described by the gluon distribution function 
f g (x). Moreover, Tq{X) coincides with Xf g (X) in the ( — > limit, and this fact suggests the 
approximation JF^(A) « Xf g (X) for small (. One should realize, however, that the electropro- 
duction amplitude is dominated by the imaginary part whose magnitude is determined by ^f(C)- 
Since the function JF^(X) rapidly varies in the region X< ( and vanishes for X = (which is not 
the case with Xf g (X)), the relation J^(C) ~ (f g (() m ay be strongly violated. This pessimistic 
expectation is not supported by a toy model for the double distribution F mod (x, y) = A(l — x—y) n , 
in which C/™ od (C) differs from J r ^ lod (() by an extra factor (1 — Q only, the latter being close to 
1 for small (. However, the structure of the double distribution F g (x,y) in general case may be 
more involved, and a detailed analysis of the interrelationship between ^(C) and C/g(C) is an 
interesting problem for future studies. 
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